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We present the first simulation of the signature of baryonic acoustic oscillations (BAO) in Lyman-Q 
forest data containing 180,000 mock quasar sight-lines. We use eight large dark-matter only simu- 
lations onto which we paint the Lyman-a field using the fluctuating Gunn-Peterson approximation. 
We argue that this approach should be sufficient for the mean signature on the scales of interest. 
Our results indicate that Lyman-o? flux provides a good tracer of the underlying dark matter field on 
large scales and that redshift space distortions are well described by a simple linear theory prescrip- 
tion. We compare Fourier and configuration space approaches to describing the signal and argue 
that configuration space statistics provide useful data compression. We also investigate the effect 
of a fiuctuating photo-ionizing background using a simplified model and find that such fluctuations 
do add smooth power on large scales. The acoustic peak position is, however, unaffected for small 
amplitude fluctuations (< 10%). Larger amplitude fluctuations make the recovery of the BAO signal 
more difficult and may degrade the achievable significance of the measurement. 



I. INTRODUCTION 



Oscillations of the baryon-photon plasma in the early 
universe, also known as Baryon Acoustic Oscillations 
(BAO), imprint a distinct signature on the clustering of 
matter P, y]- The distance that acoustic waves propa- 
gate in the first several hundred thousand years of cosmic 
evolution set a characteristic scale that is measurable as a 
distinct peak in the correlation function of matter fluctu- 
ations, or as an oscillatory pattern in the power spectrum 
of the same (see Q for a detailed description of the 
physics in modern cosmologies and for a comparison 
of Fourier and configuration space pictures). These oscil- 
lations have been traditionally measured in the Cosmic 
Microwave Background (see [f| for the latest) but with 
advent of new, large- volume galaxy redshift surveys BAO 
have been detected in galaxy clustering at low z as well 

UMMMMM- 

The use of BAO as a probe of cosmological parame- 
ters is especially enticing since the signal is at relatively 
large scale (around 150 Mpc), where the modes are still 
mostly in the linear regime. The power spectrum or cor- 
relation function can be thus computed quite accurately 
with only linear perturbation theory once one specifies 
the baryon-to-photon ratio and matter-radiation ratio. 
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which are both measured accurately from CMB acoustic 
peaks [13I . [1^ . The ability to calibrate the BAO signal 
provides a standard ruler in both the transverse and ra- 
dial directions, allowing one to measure the the angular 
diameter distance and the Hubble parameter as a func- 
tion of redshift in the clustering of matter. Absent any 
systematic errors, obtaining a high precision measure- 
ment of the distance simply requires surveying a large 
volume and locating the features in the 2-point function 
corresponding to the acoustic scale. The recent mea- 
surements of the distance scale to z ~ 0.2 — 0.4 provide 
us with complementary constraints to other large scale 
structure probes, significantly improving constraints on 
key cosmological parameters. 

The BAO technique becomes even more powerful as 
one moves to higher redshift, where the acoustic scale is 
expected to be more linear and at which more volume 
is available to be surveyed. Unfortunately tracing the 
large volumes with high fidelity, as is required by BAO 
studies, becomes increasingly expensive of telescope time 
if galaxies are used as the tracer. However, in principle 
any tracer of the mass field will do, including the neutral 
hydrogen in the inter-galactic medium (IGM) [l^, [l^ 
or galaxies 17]. Tracing neutral hydrogen in galaxies 
via its redshifted 21cm emission is a key goal for pro- 
posed future radio telescopes. However, even with cur- 
rent technology it is relatively straightforward to obtain 
a low resolution spectrum of distant quasars and study 
the Lyman-a forest of absorption lines which map the 
neutral hydrogen along the line-of-sight. At z ~ 2 — 3 
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the gas making up the IGM is in photo-ionization equi- 
hbrium, which results in a tight density-temperature re- 
lation for the absorbing material with the neutral hydro- 
gen density proportional to a power of the baryon den- 
sity [H, [3]. Since pressure forces are sub-dominant, the 
neutral hydrogen density closely traces the total matter 
density on large scales. The structure in QSO absorption 
thus traces, in a calculable way, slight fluctuations in the 
matter density of the universe back along the line-of-sight 
to the QSO, with most of the Lyman-a forest arising from 
over-densities of a few times the mean density. 

The measurability of the BAO signal in the Lyman-a 
forest has received relatively little attention in the lit- 
erature. Motivated by this and the development of the 
upcoming BOSS experiment, which will deliver unprece- 
dented number of quasar spectra probing the Lyman-a 
forest at z 2 — 3, we further develop the theory. It 
is difficult to compute the BAO signal in the forest ana- 
lytically, involving as it does a projection of a non-linear 
mapping of a non-linear density field in redshift space. 
For this reason we resort to large N-body simulations. 
Our work is an expansion of, and is complementary to, 
[20| who estimated the potential of a survey such as BOSS 
to measure the BAO scale assuming that the Lyman-a 
flux traces the dark matter on large scales, following pre- 
scriptions from [21'] . This work confirms that intuition by 
an explicit modeling of the Lyman-a forest flux, albeit in 
a simplifled manner that does not fully capture the all of 
the small-scale physics. 

The outline of the paper is as follows. We describe 
our simulations in Sec. [TTl discussing our results for the 
2-point statistics of the Lyman-a forest in Scc. lIIIBl Wc 
measure and attempt to explain the modiflcation of the 
signal by non-linear effects such as redshift-space distor- 
tions in the flux decrement correlation function. Fur- 
thermore, we make a preliminary investigation of one of 
the many systematic efl^ccts that could change the BAO 
signal ~ UV background fluctuations - in Sec. IIVI Fur- 
ther investigation is underway to thoroughly understand 
the systematic effects of the Lyman-a forest BAO signal, 
but they are beyond the scope of this paper and will be 
reported in a future publication. Finally we discuss the 
possible strategies for the detection of the BAO signal 
in the flux decrement correlation or powerspectrum and 
conclude in Sec. |Vl 



II. SIMULATIONS 

As mentioned above, it is difficult to compute the BAO 
signal in the Lyman-a forest analytically. It is also a chal- 
lenging problem numerically. Since much of the signal of 
interest comes from near mean density gas, mass resolu- 
tion is as important as force resolution. To resolve the 
Jean's scale of the gas 0(100 kpc) while simultaneously 
simulating a representative volume stretches computa- 
tional abilities, even for gravity-only simulations. We 
are thus forced to compromise. As we argue below, the 



physics governing the small-scale fluctuations should be 
approximately decoupled from the BAO scale, and inac- 
curate modeling of the small-scale physics should lead to 
smooth modifications of the flux power spectra or cor- 
relation functions which do not contain imprints of the 
acoustic scale. 

Guided by this reasoning we ran 8 particle-mesh simu- 
lations of a flat ACDM cosmology, with JIm = 0.25, JIa 
= 0.75, h = 0.75, n — 0.97 and as = 0.8. Each simulation 
evolved 3000^ particles in a 1500 /i~^Mpc box, comput- 
ing the forces on a 3000^ grid. The particle data were 
dumped at z = 2.5 and density and velocity fields gen- 
erated from the particle positions and velocities along 
a regular grid of 150^ = 22, 500 lines-of-sight using a 
spline kernel interpolation with an effective smoothing of 
250/i^^kpc. Though this smoothing is about twice the 
Jean's scale at z = 2.5, it is nearly three orders of mag- 
nitude smaller than the BAO scale providing good scale 
separation between the mis-modeled physics and the sig- 
nal of interest. 

For each line-of-sight the fluctuating Gunn-Peterson 
approximation (FGPA, [H, was used to gener- 

ate skewers of optical depth with 3, 000 pixels each. We 
assumed a temperature at mean density of 2 x 10"* K 
and an equation of state 7=1.5 Different choices 

for the slope, even an inverted equation of state, will 
quantitatively but not qualitatively change our conclu- 
sions. The optical depth included thermal broadening 
(assumed Gaussian) and skewers were generated both 
with and without peculiar velocities for the gas. The 
optical depth was scaled so that the mean transmitted 
flux F = (exp(— t)) = 0.8 [25^. For completeness, we 
also generate the skewers with dark-matter over-density 
only, so we can compare the flux statistics to those of 
the underlying mass. We work throughout with rela- 
tive fluctuations in the flux, Sp — F{x)/F — 1, so our 
fundamental data set is Sf{x) on 150^ skewers of 3,000 
pixels each. In the results presented in this paper, F 
was determined globally for each box under considera- 
tion. We have also attempted to determine F individu- 
ally for each skewer, which is closer to the observational 
situation where one has fit continuum for each individual 
quasar sight-line. As expected, the resulting two-point 
correlation functions changed by a small constant offset, 
which can be easily modeled in the analysis of the real 
data. 

In this exploratory work we neglect several higher order 
effects in order to concentrate on the underlying physics. 
First, we neglect the evolution of Lyman-a forest with 
redshift - both the structure of the IGM and the mean 
flux - and generate the skewers at z = 2.5. Second, we 
make the distant observer approximation that all skew- 
ers are parallel when in fact the comoving radial distance 
changes by about 20% across the depth of our simulation. 
This has two implications. The flrst is that the transverse 
distance between two points is to within ~ 10% of that 
obtained by assuming perfectly parallel skewers. The sec- 
ond is that in reality the line-of-sight velocities are not 
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FIG. 1: The ratio of the monopole of the flux correlation 
function in redshift space for density fields smoothed on scales 
a factor of two (filled red points) or four (blue crosses) larger 
than our fiducial case. Errors are heavily correlated and only 
meant to be indicative of the underlying uncertainties. 
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FIG. 2: The cross-correlation coefficient between the flux in 
our low and high resolution boxes, y^fj^u^hh- Filled red 
points show the result for the two low resolution boxes having 
twice the smoothing length of the high resolution box, blue 
crosses is the same for 4x smoothing length. 



perfectly parallel and that the Kaiser formalism that we 
later employ will not be exact. We will return to these 
issues in a future publication. 

There are two important effects that make our simula- 
tions inaccurate. First, the simulations do not model the 
gas, assuming it faithfully traces the dark matter on large 
scales (e.g. '2Q\). Even if we did attempt to model the gas, 
our resolution is not adequate to model the small scale 
physics; a fully converged simulation would require a res- 
olution of ~ 30 kpc, a factor of several smaller than 
we achieve [l^. Within the assumptions of the peak- 
background split the non-idealities will be absorbed into 
smooth additive and multiplicative functions and will not 
qualitatively alter our conclusions. 



To test the effect of limited resolution on the B AO scale 
we smoothed the 3D density and velocity fields of two of 
our simulations by a factor of two or four before com- 
puting the skewers. As shown in Fig. [1] the simulations 
are quite well converged, with the ratio of flux correla- 
tion functions agreeing to within 10% for a factor-of-two 
change in resolution. We expect that the missing small 
scale physics will have a larger effect on the covariance of 
the clustering, i.e. the error bars, but we do not currently 
have enough simulations to make this measurement reli- 
ably. 

We can also compute the cross-correlation between the 
flux along the same lines-of-sight in the two simulations 
with different 3D smoothings of the density and veloc- 
ity fields. This directly tests the impact of small-scale 
physics on how well the flux traces the density on large 
scales. Figure [H shows this cross correlation. We do not 
plot the errors, but the uncertainty can be gauged from 
box-to-box scatter. The correlation is nearly unity on 
large scales, indicating that the details of the small-scale 
physics do not alter the behavior of the flux on large 
scales. The signal becomes difficult to measure beyond 
100/i~^Mpc where the correlation functions are small 
and become negative. 



III. BAO SIGNAL 

A. Fourier space 

We start by discussing the signal in Fourier space, 
where we expect different k modes to be decoupled in 
the linear regime. Since the signature of BAO in Fourier 
space is a set of oscillations rather than a single peak, 
the aliasing of the modes introduced by a finite window 
function becomes severe in our problem. While this alias- 
ing is most well known from pencil-beam galaxy surveys, 
it becomes more even more acute in the Lyman-a for- 
est as the data are even more sparsely sampled. In fact, 
we were unable to obtain a measurement of the three- 
dimensional flux power spectrum from our skewers, due 
to strong mixing with the window function. The regular 
distribution of our skewers likely aggravated the prob- 
lem, but it would remain even for irregularly distributed 
skewers. 

One can avoid the aliasing issue, however, if one works 
with cross-power spectra, defined as 



qdq 



fc2 



(1) 
(2) 



where Sk is the Fourier transform of matter over-density 
or flux along a single skewer and r±_ is the transverse 
separation between skewers. 

In Figure [3] we show the cross spectrum for the matter 
overdensity and the corresponding linear predictions for 
both baryonic and baryon-less universes for a few values 
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FIG. 3: The cross power spectrum for the matter overdensity. Filled red points are measured from our simulations, the solid 
black line is the linear theory prediction, while the blue dashed line is the same for a baryon-less universe. 



of rj_ . Note that the acoustic feature not only appears as 
a series of oscillations superimposed over a smooth shape, 
but the oscillations have an increasing frequency. Meth- 
ods which involve measuring binned cross power spectra 
would need very fine binning to avoid losing signal-to- 
noise. While the full array of cross-spectra undoubtedly 
contains the acoustic information, the difficulties of ex- 
tracting it from this statistics appeared daunting, and we 
do not consider the cross power spectra any further. 



B. Configuration space 

The configuration space statistics more naturally allow 
for complex observing geometries, and irregularly sam- 
pled data while in principle having more complex covari- 
ance properties. We measured the correlation function of 
both S and Sf in both real and redshift space by brute- 
force averaging over all possible pairs of pixels in bins of 
perpendicular distance r± and parallel distance r||. 

In real space, within the context of the peak- 
background split model one can argue that the two-point 



flux statics should go as 

^ffir)^B{r)^ssir)+A{r), 



(3) 



where is the flux correlation function, ^ss is that of 
the matter and A{r) and B{r) are functions that are 
smooth on large scales. In practice we find that B(r) is 
a constant and A(r) is consistent with zero in our simu- 
lations allowing us to write 



(4) 



with b a large-scale bias. 

In redshift space, we start by modeling the large-scale 
redshift space distortions from super-cluster infall as [l^] 



E 

£=0,2.4 



(5) 



where Lf (ii) indicates the Legendre polynomial of order 
£ and [H] 

CoW = Co^(r), 

= C2 [e(r) - CM] , 



(6) 
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FIG. 4: The correlation functions, multiplied by for easier visualization, as a function of r± and r\\ in h~^Mpc. (Top) 
Correlation functions for matter over-density in real space, redshift space and the theoretical predictions for linear theory. 
(Bottom) The correlation function in real and redshift space for the flux fluctuations, Sf = F/F — 1, and (bottom right) 
the relative difference between the flux and (appropriately scaled) dark matter correlation functions in redshift space (^// — 
b^^ss)/(,f f ■ The features in this plot follow the crossing of ^// through zero, where the quantity we plot diverges. Note that 
each panel has a very different color scale. 
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FIG. 5: The monopole correlation function of the dark matter, 
multiplied by r'^, in real space (filled red points) and redshift 
space (blue crosses), the latter divided by Co- Error bars are 
estimated from the box-to-box scatter, are heavily correlated 
and only indicative of the underlying uncertainties. The linear 
theoretical prediction is plotted as the (black) solid line. 



with 

2 1 , 
Co = l + -/3+-/3^ 

3 5 

4 4 , 

C. ^ (7) 

and 

^ Jo 

f(r) = 1 rsHRis)ds. (8) 
Jo 

We use here to explicitly denote the real-space correla- 
tion function. These expressions hold for both the matter 
and the flux, with /3 ^ f = d\nS/d\na ^ fl'^^ ~ 1 in the 
former case. We leave /3 as a free parameter for the flux, 
although we expect it to be close to unity also. 

We show our main results in Figure |4l which plots 
the correlation functions (multiplied by r^) in both real 
space and redshift space, along with theoretical predic- 
tions from the simple model described above. The left 
panels show the real-space correlation function, where 
one can clearly see the expected BAO 'ridge' at r ~ 
105 /i~^Mpc. The middle panels show the correlation 
function in redshift-spacc, which is very similar to the 
left panels with the exception of a bias. The Lyman-a 
flux follows the dark matter and shares the same /3 ~ 1. 



6 




60 80 100 120 140 160 

r[/i/Mpc] 




160 




0.10 
0.08 
0.06 
0.04 
0.02 
0.00 

60 80 100 120 140 160 

r[/i/Mpc] 

FIG. 6: (Top) The flux correlation function in real (filled 
red points) and redshift (blue crosses) space. The black solid 
and and dashed lines are the scaled matter correlation func- 
tions. (Bottom) as above but divided by the equivalent mat- 
ter linear correlation function. There is no deviation from 
scale-independent bias seen in the lower panel. Errors are 
estimated from the box-to-box scatter, are heavily correlated 
and only indicative of the underlying uncertainties. 



In the upper right-hand panel we show the predictions 
calculated using the linear theory correlation function 
and super-cluster infall with /3 — 0.96, appropriate to 
the matter at z — 2.5. 

Our ability to use linear theory to describe the acous- 
tic signature at these redshifts is an important feature of 
the method. The dominant effect of non-linear cluster- 
ing is to broaden the acoustic peak, with an amplitude 
that can be estimated from the rms Zel'dovich displace- 
ment m, HO, At z = 2.5 this is about Sh'^Mpc 
in our cosmology, to be compared to the much larger in- 
trinsic width of the acoustic feature (set by the diffusion, 
or Silk, damping scale: 12 /i~^Mpc). Adding these in 
quadrature we see non-linear evolution will only change 
the peak width by 4%. Thus linear theory should accu- 
rately describe the acoustic feature in the matter at these 
redshifts. It is also interesting to note that super-cluster 
infall does not generate an elliptical contour in redshift 
space, as can be seen clearly in Eq. (jH]). 



FIG. 7: The cross-correlation coefficient between flux and 
matter over-density in real (flUed red points) and redshift 
(blue crosses) space. 



Finally, the bottom right-hand panel of Figure [3] shows 
the residuals (^// — b'^£,ss)/^ff in redshift space. The 
bright ridge is associated with the correlation function 
going though zero, which increases the fractional errors, 
but there is no large-scale radial structure in the differ- 
ence plot. 

We now focus our attention on the angularly averaged 
(or monopole) correlation function, which in the large- 
scale limit should be proportional to the real-space corre- 
lation function (Eq. [5]) . Figure [5] shows of the matter 
in real and redshift space, using only information from 
the skewers, compared to the linear theory predictions. 
We note that linear theory describes the matter very well 
on these scales and at this high redshift in both real and 
redshift space. 

Figure [S] shows the real and redshift space correlation 
functions of the flux and the ratio of the flux to matter 
correlation functions (the 'bias' squared). The implied 
bias is consistent with scale- invariant at b ~ 0.2, compa- 
rable to that obtained by [2l[. Finally, Figure[7]shows the 
cross-correlation coefficient between the flux and matter 
over-density. 



(9) 



which is unity when flux traces matter perfectly. While in 
principle r < I, when it is measured from simulations it 
can exceed unity due to noise. The simulations show that 
the variations in the flux are tracing those in the mat- 
ter remarkably well on all scales of interest, in both real 
and redshift space. This validates the idea of measuring 
the BAO feature in the Lyman-a forest, and shows that 
there is no information lost in measuring flux rather than 
matter fluctuations on these large scales even though the 
dynamic range in the flux can be drastically smaller than 
in the density field. 
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IV. FLUCTUATING PHOTO-IONIZATION 
RATE 

While our ability to measure the acoustic scale with 
galaxies has been impressively demonstrated, we still do 
not know whether this will be possible from the Lyman-a 
forest. The calculations above suggest that the signal is 
present, and the propagation of statistical errors suggest 
the measurement will be interesting [23|, however there 
are many systematic error which need to be controlled in 
order to obtain the forecast statistical precision. 

One possible cause for concern, in addition to sys- 
tematic errors in the measurements themselves, is 
non-gravitational contributions to the flux correlations. 
These could arise from hydrodynamic forces, radiative 
transfer effects, reionization heating or other departures 
from the simple FGPA assumed thus far [l^. Many of 
these effects are expected to contribute mostly on small 
scales, with no power preferentially on the acoustic scale, 
and by the arguments provided earlier should not bias 
the BAO measurement (see e.g. Figure [2]). Diagnostics 
of these non-gravitational contributions can be found in 
the higher moments of the flux [s^, [l^ and all indica- 
tions are that the forest is dominated by gravitational 
instability on large scales. However, the errors on such 
measurements are still large enough that the issue is not 
settled. 

One possible contributor to the observed power on 
large scales is fluctuations in the UV background field 
or the photo- ionization rate (F). Since the attenua- 
tion length of the IGM at z ~ 2 — 3 is large, and the 
background is thought to be dominated by rare sources 
(QSOs), r may have spatial structure on large scales. If 
the IGM is in photo-ionization equilibrium the optical 
depth, T oc F~^. Since this is the most obvious source of 
'extra' large-scale power we investigate this scenario in 
more detail. 

To begin we estimate the expected magnitude of fluc- 
tuations in the ionizing background field assuming it is 
generated by UV light from quasars [2^, [s^, [s^ . We do 
this numerically by populating each (1.5 /i~^Gpc)'^ box 
with mock quasars with luminosities following a broken 
power-law luminosity function 
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(10) 



with the parameters detailed in Table 5 of [36|. We ne- 
glect any luminosity dependent conversion from optical 
of UV luminosity in this preliminary study. Each QSO 
is assumed to emit isotropically with constant luminosity 
L, so the contribution to the photo-ionization rate from 
the ith QSO at distance r.i can taken to be 
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which neglects finite lifetimes or light-cone effects. Here 
ro is the 'attenuation length' of the IGM, which is 



TABLE I: Relative rms fluctuations in the UV background 
or photo-ionization rate, ar/{T), for different choices of the 
slopes in the QSO luminosity function and the magnitude 
limits of the QSOs we include. 



0(100 Mpc). We chose tq = lOO/i-^Mpc, close to the 
acoustic scale, to maximize the potential contamination 
effect m. 

Since tq is large, we neglect any QSO clustering in 
our model, placing the sources at random within the vol- 
ume. The spatial structure in the photo-ionization rate, 
or summed UV background, depends on the luminosity 
function, in particular on the slopes at both the faint 
and bright end. Table |T] shows some characteristic ex- 
amples, with the variations in the photo-ionization rate 
ranging from 31% to 39%. In what follows we shall take 
0.01 < L/L* < 100, a = -3.31 and /? = -1.09 as our 
fiducial model for the QSO component. The correlation 
function of F is close to constant at small scales, and falls 
dramatically beyond the attenuation length, tq. 

To the particular realization of the QSO component 
we add a uniform piece (to model the emission from faint 
AGN, galaxies and the IGM itself [11]) such that the rms 
fluctuation of the total is 5%, 10% or 20% and divide our 
original t in every pixel by the ionization rate at that 
location. The overall normalization is rescaled to set the 
mean flux, F — 0.8. 

Figure [H] shows the correlation functions of the result- 
ing fluxes. Fluctuations in F significantly affect the cor- 
relation function on large scales, increasing the level of 
power by a large factor. The additional power is however 
quite smooth, indicating that the acoustic information 
is still accessible. To see whether one can still reliably 
measure the position of the peak in such degraded data, 
we have attempted to fit the data using the following 
theoretical model 

C(r) = (^U{r) + ^G(rpeak, (Tpeak)) + I, (12) 

where 6, ft., rpcak, Cpcak and I are free parameters, <^nb is 
the linear correlation function for a baryon-less universe, 
and G(/i, cr) is a Gaussian. We do not advocate this as 
a realistic model of the correlation function, but merely 
as a convenient prescription allowing us to asses whether 
there are any biases introduced in the measurement of the 
position of the peak by the presence of UV fluctuations. 
We found best-fit point in each of the 8 boxes individ- 
ually, assuming a diagonal covariance with the plotted 
error bars, and used the scatter between these best fits 
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FIG. 8: The effect of photo-ionization rate fluctuations on the flux correlation function. The four panels show correlation 
functions for UV fluctuations of 0% (top left), 5% (top right), 10% (bottom left) and 20% (bottom right). Note that the 
vertical scale is different in each panel. The solid lines correspond to the model of Eq. 1121 



model rpeak Cpcak 

dark matter, real space 113.2 ± 1.0 16.5 ± 2.5 

flux, redshift space 113.1 ± 0.6 17.9 ± 2.9 

flux, redshfit sp., UV=5% 113.1 ± 0.6 18.6 ± 3.1 

flux, redshift sp., UV=10% 112.8 ± 0.7 20.4 ± 3.0 

flux, redshift sp., UV=20% 111.1 ± 1.1 24.4 ± 4.1 

TABLE II: Fits of the correlation function, including UV 
background fluctuations, to the model of Eg. 1121 



to asses the uncertainty in the parameters. The results 
of these fits are shown in Table |TT] and Figure [51 

These results indicate that the peak position is not af- 
fected by the presence of small UV fluctuations, but fluc- 
tuations larger than around ~ 20% start to overwhelm 
the acoustic signal and biases begin to introduced into 
the peak position recovery. It is possible that more so- 
phisticated modeling could still recover the peak, but it 
is likely that such large UV fluctuations will degrade the 
accuracy of the BAO measurement. In principle one can 
search for evidence of UV fluctuations in existing QSO 



spectra. We also find a weak evidence that the width of 
the peak broadens as UV fluctuations are increased. 



V. DISCUSSION AND CONCLUSIONS 

Baryon acoustic oscillations have become one of our 
most promising methods for determining cosmological 
distances and hence the expansion history of the Uni- 
verse. The structure in the spectrum of distant quasars, 
which is thought to trace the structure of the IGM at 
near mean density, has been suggested as a relati vely 
cheap method for measuring BAO at high redshift [isj . 
It is difficult to compute the BAO signal in the Lyman-a 
forest analytically, since it involves projection of a non- 
linear mapping of a non-linear density field in redshift 
space. Instead we have used large-volume N-body sim- 
ulations. We argued that while the small-scale physics 
in these simulations is not accurate, the overall picture 
should not change drastically with more realistic simula- 
tions. Our conclusions can be summarized as follows. 

We see clear evidence for the acoustic scale in our 
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180,000 mock spectra. While the BAO signal is present in 
both Fourier and configuration space statistics, the latter 
seem to be the better for analyzing the data. There are 
two mains reasons for this. First, the BAO feature is a 
single comparably narrow feature in the correlation func- 
tion, rather than a series of oscillations in the power spec- 
trum. Second, the complex nature by which the Lyman- 
a forest samples the underlying field makes the analysis 
considerably more subtle in Fourier space, or conversely 
the mask much easier to handle in configuration space. 

Within the peak-background split approximation, the 
Lyman-a flux follows the mass fluctuations on large 
scales. Within the FGPA the flux is a highly non-trivial 
but deterministic transformation of the underlying den- 
sity field and we find that the Lyman-a fiux correlation 
functions trace the mass with high fidelity on large scales, 
with negligible scale dependent bias. We also find that 
the cross-correlation between the flux and matter is con- 
sistent with unity. This is both surprising, given that the 
flux has limited "dynamic range" compared to the mass, 
and very encouraging. 

Given the difficulty of observing in the deep ultra- 
violet bands from the ground, Lyman-a observations are 
limited to redshifts z > 2 from ground based observa- 
tions. The upper practical limit is set by the fact that 
quasars must be bright enough to observe them with a 
survey instrument and this limits the upper redshift to 
z ^ 2 — 3. Our simulation redshift of 2 = 2.5 is thus 
at a typical value for a future instrument such as BOSS. 
Quantitative details are likely to change with redshift, 
but our main qualitative conclusions are robust, because 
there is no fundamental change in the properties of the 
intergalactic medium in the redshift range of interest. 

Similarly, our skewers were simulated at a fixed red- 
shift, while the real Lyman-a forest probes an evolving 
intergalactic medium. This will manifest itself as an ef- 
fective large-scale bias that evolves with redshift, but is 
unlikely to change abruptly. We therefore propose to 
analyse the data in discrete redshift bins that are wide 
enough to have enough Fourier modes along the line of 
sight, but are at the same time narrow enough to allow an 
approximation of a constant or linearly varying effective 
large-scale bias. 

However, a number of other systematic effects will need 
careful investigation before the BAO program can be car- 
ried out with Lyman-a forest data. We have begun by 
investigating the effect of spatial structure in the photo- 
ionization field, which modulates the optical depth on 
scales comparable to the acoustic feature. We find that 
such a modulation does not introduce any confusion into 
the determination of the acoustic scale, though it does 
modulate the total amplitude of the flux correlation func- 
tion. This is not too surprising, the photo-ionization rate 
does not know to add power at 105/i~^Mpc differently 
than at 104 or 106, but is also encouraging. 



As the simulations become increasingly realistic we 
may begin to see some scale-dependence in the bias. 
We therefore suggest that one should model BAO in the 
Lyman-a forest, by measuring a set of bands in the cor- 
relation function, together with a linear /? parameter. 
This could for example be achieved using an (optimal) 
quadratic estimator. The resulting /3 should be close to 
that of the matter, i.e. close to unity. In order to be 
conservative, the resulting correlation function could be 
compared with the theoretical predictions for the matter 
correlation function using the formalism in Equation [3] 
with A(r) and B(r) being smooth, one or two parameter 
functions. Sufficiently accommodating functions can also 
protect against unexpected systematics. 

Significantly more work needs to be done before this 
program can be executed and there are several important 
effects that need to be understood. First, quasar contin- 
uum fiuctuations can potentially give rise to large-scale 
fiuctuations that might have a preferred scale. The fact 
that we are cross-correlating different spectra alleviates 
the problem to some extent, but may not solve it. Simi- 
larly, metal contamination can also give rise to spurious 
correlations. If fluctuations in the UV background, or 
photo-ionization rate, are large they may imprint large- 
scale power in the flux correlation function which can 
overwhelm the acoustic signal and reduce the sensitivity 
of the measurement. 

Many of these issues can be dealt with using two ba- 
sic approaches. The first is to model them. For exam- 
ple, pixel pairs which are thought to be contaminated 
by a metal doublet could be 'blinded' within an opti- 
mal quadratic estimator. A second method is to split 
the full sample into all possible subsamples. For exam- 
ple, quasars of the same absolute magnitude or the same 
source redshift are more likely to have similar continuum 
fiuctuations and one could check the results by splitting 
the sample into subsamples according to QSO magnitude 
or redshift. Our results provide added impetus to further 
develop these promising ideas. 
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